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Abstract 

We investigate simplified models of computer data networks and examine how the 
introduction of additional random links influences the performance of these networks. 
In general, the impact of additional random links on the performance of the network 
strongly depends on the routing algorithm used in the network. Significant performance 
gains can be achieved if the routing is based on "geometrical distance" or shortest path 
reduced table routing. With shortest path full table routing degradation of performance 
is observed. 
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1. Introduction 

Models of computer data networks have attracted much attention in recent years. Generally, 
these models assume either regular topology of the network, in the form of a square lattice 
]!], [2|, |7J or a binary Cayley tree [ II], or random graph topology ||. On the other hand, 



it has been recently demonstrated that many technological, biological, and social networks 
are neither completely regular nor completely random, being somewhere between these two 
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extremes ||. Regular lattices rewired to introduce small amount of random connections, 
termed "small world networks" , offer many advantages over purely regular or purely random 
topologies. In particular, models of dynamical systems based on "small world" lattices can 
often exhibit enhanced signal propagation capabilities, as observed in epidemic models or 
multi-player prisoners dilemma games played on such lattices [|J. 

The purpose of this work is to investigate if a similar effect can be achieved in simple 
data network models, although our approach is different than the approach taken in [[|. 
We do not rewire the network, but rather examine how introduction of additional random 
links influences its performance. This is motivated by the question whether introduction of 
additional links can help to decongest an existing network. 

2. Network Models Definitions 

The purpose of the network is to transmit messages from points of origin to destination 
points. In our model, we will assume that the entire message is contained in a single "capsule" 
of information, which, by analogy to packet-switching networks, will be simply called a 
packet. In a real packet-switching network, a single packet carries the information "payload" , 
and some additional information related to the internal structure of the network. We will 
ignore the information "payload" entirely, and assume that the packet carries only two pieces 
of information: time of its creation and the destination address. 

Our simulated network consists of a number of interconnected nodes. Each node can 
perform two functions: of a hosts, meaning that it can generate and receive messages, and 
of a router (message processor), meaning that it can store and forward messages. Packets 
are created and moved according to a discrete time parallel algorithm. The structure of the 
considered networks and the update algorithm will be described in subsections which follow. 

2.1. Connection Topologies 

We will consider two types of network connection topologies: a two-dimensional square lattice 
C = C(L) and a two-dimensional square lattice £ with additional links added randomly, 
denoted by Ci = £i(L). The value of the subscript / gives a number of an extra links in a 
network and the value of L gives a number of nodes in the horizontal and vertical direction 
of the lattice C. The lattices £ and L\ with periodic boundary conditions will be denoted 
by C p and £f, respectively, and with non-periodic boundary conditions by C np and £™ p , 
respectively. Hence, with this notation, Cq = C p and Cq P = C np . Most of our simulations 
will be performed on lattices with periodic boundary conditions. 

The network hosts and routers are located at nodes of the lattice L. The position of each 
node on a lattice £ is described by a discrete space variable r, such that 

r = ic x +jc y , (1) 

where c x , c y are Cartesian unit vectors, and i, j = 1, . . . , L. For each node r we denoted by 
C(r) the set of all nodes directly connected with the node r. Hence, for each r G C p , the set 
C(r) is of the form 

C(r) = {r - c x ,r + c x ,r - c y ,r + c y }. (2) 
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In this case, the node r is connected with its four nearest neighbours. However, for lattices 
with non-periodic boundary conditions C np or square lattices with additional links added 
randomly C p and the form of the set C(r) can be different, for some nodes r, from the 
one in (0). For example, for nodes r on the boundary of a lattice C np the set C(r) can contain 
two or three elements only, depending on where a node r is located on the boundary. In 
the case of a lattice C p or £™ p the set C(r) can contain many non-nearest neighbours nodes 
depending on a number of additional links which originate from the node r. 

The extra links are constructed using the following procedure. We first select randomly 
a node ri on a square lattice C p or C np . Next, we select randomly another node r 2 , different 
from the node ri, and connect these two nodes with direct communication link. By repeat- 
ing this procedure independently I times we obtain a lattice C p or respectively, with 
additional / random links. It can happen that the nodes ri and r 2 can be selected again to 
form a new link. Hence, in the network there can be several links connecting directly the 
same nodes. We want to emphasize that all the connections in our models are static, during 
the simulation period they do not change. Additional random links are added before the 
simulation starts, and remain unchanged. 

In the networks considered here, each node maintains a queue of unlimited length where 
the arriving packets are stored. The number of packets in the queue at a node r at time step 
k will be denoted by n(r, k), while the total number of packets in the system at time step k 
will be denoted by N(k), 

N(k)=Y,n(r,k). (3) 

Packets stored in queues, at individual lattice nodes, must be delivered to their destination 
addresses. To assess how far a given packet is from its destination, we introduce the concept 
of distance between nodes. Depending on a network connection topology we will use three 
metric functions to compute the distance between nodes r x = («i,ji) and r 2 = (^j^)- 
Namely, we will use 

1. for lattices with non-periodic boundary conditions "Manhattan" metric 

d M (ri,r 2 ) = \i 2 - h\ + \j 2 ~ (4) 



2. for lattices with periodic boundary conditions periodic "Manhattan" metric 



d PM (r 1 ,r 2 ) = L 



\i2 ~ ill 



\h - jil 



(5) 



3. and regardless of boundary conditions the "shortest path" metric <isp(ri, r 2 ) defined 
as the number of links in the shortest path joining ri and r 2 . By "shortest" we mean 
the path with the smallest number of links. 

Metric dpu will be used on lattices C p , while du will be used on lattices C np . Note that 
on a square lattice with no extra links dsp(Fi, r 2 ) = d M (r 1 , r 2 ), and dgpfti, r 2 ) = d PM {vi, r 2 ) 
for lattices with non-periodic and periodic boundary conditions, respectively. Furthermore, 
for each r G Cf, where k = p or np and I G {0, 1, ...} 
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C(r)={x6^ SP (x,r) = l}. (6) 

When irregularities such as extra links are present, dsp(vi, r 2 ) can be computed using one of 
the well known algorithms. In our simulations, we used shortest path backward tree algorithm 



2.2. Routing Algorithms 

The dynamics of the networks are governed by the parallel update algorithms shown Figure 
0], similar to the algorithm used in 0. We start with an empty queue at each node, and 
with discrete time clock k set to zero. Then, the following actions are performed in sequel: 

1. At each node, independently of the others, a packet is created with probability A. Its 
destination address is randomly selected among all other nodes in the network with 
uniform probability distribution. The newly created packet is placed at the end of the 
queue. 

2. At each node, one packet (or none, if the local queue is empty) is picked up from the 
top of the queue and forwarded to one of its neighboring sites according to a one of 
the routing algorithms to be described below. Upon arrival, the packet is placed at 
the end of the appropriate queue. If several packets arrive to a given node at the same 
time, then they are placed at the end of the queue in a random order. When a packet 
arrives to its destination node, it is immediately destroyed. 

3. k is incremented by 1. 

This sequence of events, which constitutes a single time step update, is then repeated 
arbitrary number of times. The state of the network is observed after sub-step 3 (clock 
increase), but before sub-step 1 (creation of new packets). In order to explain the routing 
algorithms mentioned in sub-step 2, we will first describe one of its simplified versions. 

Let us assume that we measure distance using some metric d, where d could be any of the 
previously defined metrics du , dpM, or dsp. To decide where to forward a packet located at 
a node r with the destination address r<j, two steps are performed: 

1. From sites directly connected to r, we select sites which are closest to the destination 

of the packet. More formally, we construct a set A^lr) such that 

Aoo(r) = {a G C(r) : d(a.,r d ) = min d(x, r d )} (7) 

xeC(r) 

2. From A^r), we select a site which has the smallest queue size. If there are several 
such sites, then we select one of them randomly with uniform probability distribution. 
The packet is forwarded to this site. Using a formal notation again, we could say that 
the packet is forwarded to a site selected randomly and uniformly from elements of a 
set -Boo(r) defined as 

-Boo(r) = {a G A^r) : n(a, k) = min n(x, k)}. (8) 
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To summarize, the routing algorithm R^ described above sends the packet to a site 
which is closest to the destination (in the sense of the metric d) , and if there are several such 
sites, then it selects from them the one with the smallest queue. If there is still more than 
one such node, random selection takes place. It is clear that each packet routed according 
to the algorithm Roo will travel to its destination taking the shortest possible path (shortest 
in the sense of the metric d, not necessarily in terms of the number of time steps required 
to reach the destination). In real networks, this does not always happen. In order to allow 
packets to take alternative routes, not necessarily shortest path routes, we will introduce a 
small modification to the routing algorithm Roo described above. 

The modified algorithm R m , for each node r, will use instead of the set A^t) a set A m (r) 
defined as follows. In the construction of the set A m (r) instead of minimizing distance to 
the destination d(x, r d ), as it was done in (|7|), we will minimize m (cf(x, r^)), where 

® m ^) | m, otherwise, ^ 
for a given integer m. Thus, the definition of the set A m (r) is 

A m {r) = {ae C(r) : 6 m (d(a, r d )) = min m (d(x, r d ))} (10) 

xGC(r) 

The above modification is equivalent to saying that nodes which are further than m distance 
units from the destination are treated by the routing algorithm as if they were exactly m 
units away from the destination. If a packet is at a node r such that all nodes directly linked 
with r are further than m units from its destination, then the packet will be forwarded to a 
site selected randomly and uniformly from the subset of C (r) containing the nodes with the 
smallest queue size in the set C(r). It can happen that the selected site can be further away 
from the destination than the node r. 

Therefore, introduction of the cutoff parameter m adds more randomness to the network 
dynamics. One could also say that the destination attracts packets, but this attractive 
interaction has a finite range m: packets further away than m units from the destination are 
not being attracted. 

It is also possible to relate various values of the cutoff parameter m to different types of 
routing schemes used in real packet-switching networks. Assume that each node r maintains 
a table containing all possible values of d(x, r<j), for all possible destinations and all 
nodes x e C(r), and that packets are routed according to this table by selecting nodes 
minimizing distance, measured in the metric d, travelled by a packet from its origin to its 
destination. Such a routing scheme is called table-driven routing || and it is equivalent to 
the routing algorithm R^. In this case, construction of the set A»( r ) would require looking 
up appropriate entries in the stored table. 

Let us now define D max to be the largest possible distance between two nodes in the 
network. When m < D max , then for a given x, we need to store values of d(x, r^) only for 
nodes which are less than m units of distance away - for all other nodes distance does not 
matter, since it will be treated as m by the routing algorithm. Hence, at each node r the 
routing table to be stored is smaller than in the case when m = D max . The routing scheme 
based on this smaller routing table is called the reduced table routing algorithm || and it is 
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equivalent to the routing algorithm R m . In the case when m = D max the routing algorithm 
R m Roo- 

Finally, when m — 1, the distances between hosts and destinations are not considered 
in the routing process of packets. Therefore, there is no need to store any table of possible 
paths at nodes of the network. This case corresponds to the table-free routing algorithm || 
in which packets are routed randomly. Hence, this algorithm can send packets on circuitous 
and long routes to their destinations. The analysis of this routing algorithm has been done 
in 0, where some analytical results are also presented. At present such results are not 
available for routing algorithms with m > 1. 



3. Performance of networks with square lattice connection 
topology (£ p ,d PM ), {Cl P ,d M ) 

In order to asses the performance of a network, graphs of delay as a function of presented 
load are frequently used in network performance literature ||. In our case, delay r will be 
defined as the number of time steps elapsed from the creation of a packet to its delivery to 
the destination address. We will also use average delay r(k), where the average is taken over 
all packets delivered to their destination from the beginning of the simulation (k = 0) up to 
time k. Probability of a packet creation A will be used as a measure of a presented load. 



3.1. Full table routing 

We assume that the network topology is a square lattice £q(L) with dpM metric and the 
network routing algorithm is the full table routing algorithm R m , with m = D max , i.e. 
R m = Roc. Figure Qa shows graphs of the average delay r(k) versus presented load A, as 
measured during simulation performed on a lattice £q(50). The three curves shown there 
correspond to different times. It is clear that beyond a certain critical value of A = A c , 
the average delay drastically increases. Moreover, the average delay grows with time, which 
suggests that for A > A c there is no equilibrium state. In fact, when A > A c , a typical queue 
size and consequently, the number of packets in the system N(k), grows without bounds, as 
shown in Figure and Figure [| 

It is possible to find an approximate value of the critical load A c by the following argument. 
For A < A c , the system reaches steady state, and in the steady state the number of packets 
created per unit time (given by L 2 X) must be equal to the number of packets delivered per 
unit time. Since the average time spent in the system by a packet is r(k), we can reasonably 
assume that N(k)/r(k) packets are delivered to their destinations per unit time, hence 

if- 2 - 

This relationship, known as Little's law in queuing theory ||, holds only below the critical 
point, as shown in Figure 

For the routing algorithm R m , with m = D max , when the number of packets in the 
network is small, an individual packet is always routed in such a way that it follows the 
shortest path to its destination avoiding all occupied nodes. This means that for small 
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N(k), the average packet delay is approximately equal to average distance from the packet's 
origin to its destination, which will be called "free packet" delay Tq 



to 



(12) 



ri,r 2 



After some algebra, this leads to 



1 L_1 



U,i2,jl,j2=0 
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(13) 



Obviously, when the load increases, at some point the number of packets in the network will 
be so large that it would not be possible to find a route to a destination completely avoiding 
other packets. Assuming that packets are approximately uniformly distributed over the 
entire lattice, this will happen when all sites are occupied, i.e. when N(k) = L 2 . Using (|TT|) 
this gives an estimate of A c : 

A, = =r- (14) 

TO 



For £q(50) we obtain r = 25 and A c = 0.04, in good agreement with the value obtained 
from simulations A c = 0.039 ± 0.001. 

Quite similar calculations can be performed for a square lattice with non-periodic bound- 
ary £q P (L). In this case, 



f5> 



M(ri,r 2 J 



ri,r 2 



2L 2 -1 
3 L 



(15) 



yielding A c = 0.03 for £7(50). The measured value of A c for £q P (50) is 0.020 ± 0.001, i.e., 
much lower. The discrepancy is mainly due to the fact that for the lattices Cq P (L) packets 
are not uniformly distributed on the lattice, having a tendency to cluster at the center. 
Consequently, jamming occurs earlier than one would expect assuming uniform distribution 
of packets. 



3.2. Partial table routing 

Decrease in value of the cutoff parameter m has a profound effect on the critical load. Smaller 
m means that packets which are located further than m links from their destination move 
with a high degree of randomness, and as a result, their average delay is larger. This increase 
of the delay can be also seen in a plot of a single packet delay as a function of m (Figure |5|). 
While values of m close to D max do not significantly change To, values of m close to 1 result 
in an increase of r by up to two orders of magnitude. 



4. Performance of networks based on square lattices [C p t , cIpm), 
(jC^ p , (1m) with additional / random links 
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4.1. Full table routing 

Let us now consider the network dynamics governed by the routing algorithm taking 
place on lattices £f and £™ p which in addition to normal nearest neighbor connections, feature 
I additional links, where I > 0. Figure |6| shows how addition of random links changes the 
graph of delay vs. load for both non-periodic and periodic case. We are still using cLm and 
dpM metric for Cf p and £f lattice, respectively, which means that the distance between two 
points i - ! and r 2 is still computed using d M (rx,r 2 ) or <ip A f(r 1 ,r 2 ) metric, respectively, even 
if ri and r 2 are directly connected by some extra link. 

As expected, addition of extra links improves performance of the network, shifting the 
critical point A c to the right (see Figures ^| and 0). This means that the network can carry 
more load without experiencing congestion. Performance improvement is more pronounced 
for lattices with non-periodic boundaries, as shown in Figure |f[ For example, by adding 100 
random links to 50 x 50 lattice, which increases total number of links by 2%, we increase the 
critical load by over 25%. Increasing the number of links by 8% doubles the critical load. 
This can be attributed to the fact that some packets can bypass congested central area by 
using "shortcuts" , and their delay decreases not only because they have shorter distance to 
travel, but also because they have avoided congestion. In the case of lattices with periodic 
boundaries, packets are more uniformly spread even in the presence of additional random 
links. Thus, the performance improvement is only caused by the decrease in the distance 
traveled, but not by bypassing congestion, since congestions are also uniformly spread in the 
case of lattices with periodic boundaries. 

In the remainder of this article, we will focus our discussion on lattices with periodic 
boundaries only. 



4.2. Partial table routing 

When a routing algorithm R m , with m < D max , is used, additional random links can signif- 
icantly increase critical load, and the relative performance gain is much larger than in the 
case of the full table routing algorithm. Figure ^]a,b shows the relative change of the critical 
load, defined by 

AA C _ \ c (m,l) - A c (m,0) 

A c ~ A c (m,0) ' 1 j 

where A c (m, I) denotes the critical load at a given m and I, for two different values of m, 
m = 50 and m = 20. One can immediately notice that the impact of additional links on 
performance of the network is much stronger in the case of partial table routing (m = 20) 
than in the case of full table routing (m = 50). For example, about 50 extra links are 
sufficient to double the critical load corresponding to m = 20, while the same number of 
links has almost negligible impact on the critical load when m = 50. 
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5. Performance of networks based on square lattice (C^dsp) with 
additional / random links and dsp metric 

As stated before, for a square lattice without additional links, metric d$p is identical to 
c?m or dpM metric. This is no longer true for a square lattice with additional random links. 
Routing based on dsp metric fully utilizes shortcuts provided by additional links, significantly 
decreasing "free packet" delay r . One would expect that a decrease in "free packet" delay 
will decrease also average delay, as it was in the case of networks with du and dpM metric. 
In reality, we observe just opposite effect (Figure ||). 

5.1. Full table routing 

Figure [8]c shows how the critical load A c (m, I) changes when additional random links are 
introduced. This is shown for the network dynamics governed by Roo, i.e. the full table 
routing algorithm with m = D max , on a square lattice with periodic boundaries One can 
clearly see that if the number of additional random links I is below some critical value 
l c (m) the critical load X c (D max , I) is actually smaller than \ c (D max , 0), in spite of increased 
connectivity between nodes of the network. The performance of the network is at its worst 
when just a few additional random links are added. However, it improves with the increase 
of a number I of additional random links and at some critical value l c (m) it becomes the same 
as of the network without any additional random link. When the number I of additional 
random links is greater than l c (m) an improvement in the network performance is observed. 
For the network £5 00 (50) the critical load A c (50, 500) is almost equal to the critical load 
A c (50,0) of the network £q(50), i.e. A c (50, 500) ~ A c (50,0), and the improvement of the 
£f (50) network performance is observed for I greater than 500. 

This rather unexpected phenomenon can be understood as follows. When additional 
links are introduced, and their number is less than l c {m), they provide a shortcut between 
distant parts of the network. Since packets are forwarded to their destinations via the 
shortest path, it often happens that one link serves as a shortcut for many packets from the 
neighborhood. One could say that additional links "attract" most of the traffic and quickly 
become congested, even though sites which are not close to extra links are almost empty. 
This is well illustrated in Figure |], which shows snapshots of dynamics of the network with 
-Roo routing algorithm, 50 x 50 nodes and periodic boundary conditions. The presented load 
is A(50, 0) = 0.025, just below the critical value A c (50, 0) = 0.028. If there are no additional 
random links the network dynamics remains in the steady state, as is showed by the left 
column of Figure |9|. The number of packets in the network fluctuates slightly over time, 
but remains at the same level: at k — 100 there are 2106 packets in the network, while at 
k = 1000 there are 2127 packets. When additional 100 random links are introduced, which 
is less than Z c (50), keeping all other network parameters unchanged, the network dynamics 
enters the congested phase, as it is illustrated by the right column of Figure [| The number 
of packets in the network increases rapidly over time from 2238 packets in the network at 
k = 100, to 14990 packets at k = 1000. Congestions occur mainly at inputs and exists from 
the extra links. At these nodes the queue sizes are substantially larger than in other nodes 
of the network. This is illustrated by the dark spots in the figures of the right column of 
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Figure |[ 



5.2. Partial table routing 

The performance of a network changes from the one described above when the value of the 
cutoff parameter m is less than D max , i.e. m < D max . For the routing algorithm i? 20 applied 
on the lattice £f (50) for various values of the parameter / the performance of the network is 
shown on Figure §d. From this figure we observe that the critical value l c (m) below which the 
performance of the network with the routing algorithm i? 20 is worse than the performance 
of the network without random links added is rather low. This value is lower than the 
corresponding value when Roc routing algorithm has been used. For example, adding more 
than about 15 links increases the critical load A c . Adding about 50 links, just a 2% increment 
in the number of links, increases the critical load by 100% ! The performance of the network 
improves significantly further with the increase of the number of random links. 

The explanation of this behavior is straightforward. As we have already mentioned, for 
a regular square lattice without random links added, the values of the critical load A c (m, 0) 
strongly decrease with decrease of the cutoff parameter m. This results from the fact that 
packets which are further away than m units from their destinations are not being attracted 
to the destinations and travel randomly through the network. However, addition of random 
links significantly decreases the average distances between network nodes in the metric dsp- 
Therefore, it does not matter what is the exact value of the cutoff parameter m. Most 
of the time distances in the metric dsp are way below m and packets can be attracted to 
their destinations much faster. This attraction increases with the increase in the number of 
additional random links. Hence, when additional links are present, the critical load X c (m, I) 
is not very much dependent on m, unless m is very small and increases with increase in value 
of the number / of the random links added. 

Let A c (m, /) be the critical load of a network with the cutoff parameter m and with / 
extra links added. Let rrii = D max and let m 2 be smaller than mi, but not too close to 1, 
for example, m x = 50, m 2 = 20, as in Figure The performance of the networks with the 
routing algorithms R mi and R m2 , as in the Figure || can be summarized as follows 

l c (m 2 ) < i c (mi), 
A c (m 2 ,0) < A c (mi,0), 
\ c {m h l) < A c (m;,0), 

for i — 1,2 when I < l c (mi) and 

A c (m,j, 0) < \ c (rrii,l), 
for i = 1, 2 when I > ! c (mj), and for sufficiently large / 

A c (mi,/) « A c (m 2 ,Z). 

Figure shows how the introduction of 100 additional links to a network with 50 x 50 
nodes, routing algorithm i? 20 and presented load A = 0.008 (which is above the critical load 
for / = 0), affects the network dynamics. When / = , the number of packets in the network 
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grows with time, from 1248 at k = 100 to 6750 at k = 1000, indicating that the system is in 
the congested state. The left column of Figure |10| shows that the queue sizes grow almost 
uniformly over all nodes of the network. Hence, the congestion is distributed uniformly over 
all nodes of the network. However, when I = 100 additional random links are introduced, 
the right column of Figure [H] shows that congestion is eliminated. The number of packets in 
the network remains almost steady and it fluctuates around 230 (Figure [10]). If an occasional 
small congestion occurs near the entrance to one of the shortcuts, it quickly disappears. For 
example, dark square visible in the right column of Figure W at k — 100 is not visible at 
k = 1000, demonstrating that local congestions are not permanent. 



6. Conclusion 



We found that the impact of additional random links on the performance of the network 
strongly depends on the routing scheme used in the network. Critical load of a network can be 
notably improved if the routing is based on a "geometrical distance" . Adding small number 
of additional links can decrease the average delay and shift the transition to the jammed 
phase toward higher load values. This, in general, is not true for routing schemes based 
on the "shortest path metric". In this case, if the number of additional links is small, one 
can actually observe degradation of performance: many packets attempt to utilize shortcuts 
introduced by additional links, causing congestion which in effect pushes the network to a 
jammed phase. Reduced table routing can, to some extent, eliminate this problem. If packets 
located further than m links from the destination are routed randomly (other packets taking 
the shortest possible path), performance gains obtained by adding sufficient number of extra 
links can be quite significant. 

In order to relate our findings to data network protocols used in practice, more research 
is clearly needed. In particular, congestion control mechanisms built into protocols such as 
TCP/IP will certainly affect phenomena reported here. This issue, as well as other possi- 
ble modifications of the model, is currently under investigation. Furthermore, the authors 
believe that some of the issues raised in [fTO], related to self-similar traffic modeling and 
analysis, and performance modeling of modern high-speed networks can be addressed by the 
methodology of this paper. 
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START 



c 



Procedure CREATE 



Generate random number 
qe[0,l) 




Randomly select one 
lattice node r H 



Create packet with 
destination address r d , 
place it at the end of Q(r) 



c 



RETURN 



1 




k:=0 

For every r, 
Q(r)=0 (empty queue) 






simultaneously apply 
procedure CREATE 
to all nodes r 


> 


> 


simultaneously apply 
procedure ROUTE 
to all nodes r 


' 




k:= 


k+1 




Pickup one packet from 

the top of queue Q(r), 
determine its addrerss r,. 



Construct set C(r) of all 
sites directly connected to 
r 



Find A m (r), set of sites 
belonging to C(r) which 
are closes to r H 



Construct set B m (r) 
containing sites from A m (r) 
with smallest queue size 



Randomly select one 
element in B m (r) 
and denote it a 



No 



Place the packet 
at the end of queue Q(a) 




Destroy the packet 
(packet delivered 
to its destination) 



RETURN 



Figure 1: Network update algorithm. Symbol Q(r) denotes the queue at node r. 
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Figure 2: (a) Average lifetime of a packet r(k) as a function of A, for £q(50) with dpM and 
m = D max after k = 1000 (o), k — 1500 (•), and k = 2000 (*) iterations, (b) Number of 
packets in the system N(k) as a function of A after k = 1000 (o), k — 1500 (•), and k = 2000 
(*) iterations. 
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Figure 3: Number of packets in the system N(k) for subcritical and supercritical values of 
A (A = 0.035 and A = 0.042, respectively). £q(50) with dpu metric and m = D max . 
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Figure 4: Verification of Little's law for a lattice £q(50) with dpM and m 
k = 1500. Continuous line corresponds to L 2 ) = A. 
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Figure 5: Free packet delay t as a function of m for a lattice Cq(50) with dpM metric. 
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Figure 6: (a) Average lifetime of a packet r(k) as a function of A for k = 1500 and for 
£™ p (50) lattice with du metric and m = D max , and with the number of extra random links 
I = (o), I = 100 (•), I = 200 (*), and I = 400 (+). (b) The same plot for £f(50) lattice 
with dpM metric and m = D max . 
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Figure 7: Critical load A c as a function of a number of extra links / for the lattice 50 x 50 
with periodic (o) and non-periodic (•) boundaries, using dpu and d M metric, respectively. 
In both cases, m = D max . 
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(b) dpM metric, m = 20 



Li i i i 



1000 




1000 



Figure 8: Change in critical load for a lattice £q(50) using metrics dpu and d SP with m = 20 
and m = 50. Vertical axis corresponds to (A c (m, I) — X c (m, 0))/A c (m, 0), where A c (m, I) is a 
critical load at a given m and 
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Figure 9: Comparison of dynamics of the network £f(50) with dsp and Rqo algorithm for 
I = (left column) and I = 100 (right column). Queue sizes are represented as shades of 
gray, from the highest queue size of 20 or more represented by black color to the empty 
queue represented by white color. In order to preserve clarity, additional links are shown for 
k — only. 
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Figure 10: Comparison of dynamics of the network C\ (50) with rf^p and R 2 o algorithm for 
/ = (left column) and I = 100 (right column). Queue sizes are represented as shades of 
gray, from the highest queue size of 20 or more represented by black color to the empty 
queue represented by white color. In order to preserve clarity, additional links are shown for 
k — only. 
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